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A three-parameter equation of state (EOS) without physically incorrect oscillations is proposed based on the 
generalized Lennard-Jones (GLJ) potential and the approach in developing linear isotherm regularity (LIR) EOS 
of Parsafarand Mason U- Phys. Chem., 1994, 49, 3049]. The proposed (GLIR) EOS can include the LIR EOS therein 
as a special case. The three-parameter GLIR, Parsafar and Mason (PM) [Phys. Rev. B, 1994, 49, 3049], Shanker, 
Singh and Kushwah (SSK) [Physica B, 1997, 229, 419], Parsafar, Spohr and Patey (PSP) []. Phys. Chem. B, 2009, 
113, 11980], and reformulated PM and SSK EOSs are applied to 30 metallic solids within wide pressure ranges. 
It is shown that the PM, PMR and PSP EOSs for most solids, and the SSK and SSKR EOSs for several solids, have 
physically incorrect turning points, and pressure becomes negative at high enough pressure. The GLIR EOS is 
capable not only of overcoming the problem existing in other five EOSs where the pressure becomes negative 
at high pressure, but also gives results superior to other EOSs. 
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1. Introduction 

The equation of state (EOS) describes the relationships of a system among thermodynamic vari- 
ables such as pressure, temperature and volume, which plays an important role in many fields, such 
as condensed-matter physics and geophysics. The JVIurnaghan |1] and Birch |2, 3] EOSs for solids are 
widely used in geophysics. Since Rose et al. |4] proposed in 1986 that there exists a universal EOS (UEOS) 
being valid for all types of sohds through analyzing the energy band data, a lot of forms of UEOS have 
been put forward I5|-[l6ll. However, some of them have been applied to the study of thermodynamic prop- 
erties of liquids | 6-19], while the traditional Tait EOS has also been used as universal equation both for 
solids 111,12] and liquids l lillij. 

In 1994, Parsafar and Mason (PM) proposed the following EOS by using a series expansion of internal 
energy (nl] 



Here, Vo is the volume at zero pressure. Q, Q, C2 are three coefficients in the PM EOS. In 1997, Shanker, 
Singh and Kushwah (SSK) proposed the following EOS fl2l[lil 

P^Do + Dif-^l+Daf-^l , (2) 



where Dq, Di, D2 are three coefficients in the SSK EOS. It can be seen that the SSK EOS can be expressed 
as volume-anal5?tic and pressure-analytic forms. 
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In 1994, Parsafar and Mason proposed the following linear isotherm regularity (LIR) EOS for gases 
and liquids based on the Lennard-Jones (LJ) (12-6) potential II jl 



(Z-l)y ^A,.M[-j . (3) 

Here, Z is the compressibility factor, which is equal to PVIRT. The upper density limit of LIR 1 14] is less 
certain but seems to be the freezing line for liquids (T < Tc) and at least about twice the Boyle density 
for supercritical fluids. LIR EOS has been extended to mixtures |15] and to other forms llfrtlSll through 
adopting different potential functions, including the exponential-6 |16], LJ (6-3) I17j, LJ (885-4) and 
LJ (12-6-3) 119] potentials. 

Recently, Parsafar, Spohr and Patey (PSP) jlill . extended the equation (3) to the following form with 
three parameters based on an effective near-neighbor pair interaction of an LJ (12-6-3) potential 



The PSP EOS can be equivalently reformulated as truncated Virial form 

RT Qi Q2 Q3 
V 

Parsafar et al. fli*] claimed that the PSP EOS l|4) can be applied to aU fluids and solids, and their application 
for solids [19] does not reveal any pressure or temperature limitations. 

However, we noticed that the PM EOS (T) and the PSP EOS (|4), l[5) are physically wrong at high pressure 
conditions for some solids. This is because the coefficients C2 in equation (1) and A2 in equation (4) should 
be positive for all solids to ensure a physically correct tendency at high pressure, P ^ 00 as 1/ ^ 0. 
However, the values of C2 for most solids studied in this paper are negative; and the values of A2 for 
solids NaCl and CaO studied by Parsafar et al. |19], and for most solids studied in this paper are also 
negative. This leads to an unphysical tendency, P -00 as y ^ 0. 

The incorrect tendency makes the PM Jill] arid PSP E3l EOSs inapplicable to high pressure conditions. 
We may preliminarily analyze the reason for the failure of two EOSs as follows. Holzapfel 1 20] has pointed 
out that the limitation of an EOS as the volume tends to zero, should be the Tomas-Fermi (TF) model, P oc 
V~^'^. The repulsion terms in PM ^ and PSP (H EOSs are, P oc and P oc V~^, respectively. Their 
exponent numbers 4 and 5 are far larger than 5/3, and are too hard for solids. In order to fit experimental 
P-V data at low and middle pressure ranges, the optimized C2 and A3 should take on negative values. 

In this work, we propose generalized LIR (GLIR) EOS based on a near-neighbor pair potential of the 
extended Lennard-Jones (mi,«i) type. The GLIR contains three parameters and can overcome the defect 
appearing in the PM EOS Q) and PSP EOS (5). In section 2, the three-parameter GLIR EOS is proposed. 
In section 3, equations l[l) and (2) and their modified version, PSP EOS (4) and the GLIR EOS are applied 
to twenty solids within wide pressure ranges of hundreds GPa and at ambient temperature, the results 
being analyzed and discussed. In section 4, the conclusion is presented. 



2. Analytic equations of state 

We adopt the effective pair interaction of an extended Lennard-Jones {mi,ni) type potential flQlEvi 

(6) 

mi - Ml L \ r I \ r I i 

It is well known that the effective potentials for metals usually have oscillating tails due to Friedel oscil- 
lations of electron density, and the Lennard-Jones (LJ) potentials are not really appropriate for correct 
reproduction of the energetics of metals. However, many works |10-19] have shown that the LJ poten- 
tials can mimic many properties of metals in some compression ranges. By adopting the nearest neighbor 
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assumption the total configurational energy of a solid is 



NSeo [ /re A mi 



2 (mi 

NSeo 
2{m\ - «i) 



mi 



mi/3 



mi 



ni/3 



(7) 



where V - o'lj, Vb = (re)^ Ij, a is the nearest neighbor distance, and 8 is the mean coordination num- 
ber |'lO,'l^. Following Parsafar and Mason (8||, the internal pressure can be obtained by the derivative of 
equation (6) 



dlJ_ 



miniNSeo 



mi/3+l 



6(mi - ni) 

Let us substitute the equation Q into the following internal energy equation 

{-] ■ 

[dVlr 

After integration, we derive the equation 



P= — + ^1 — 
V \V 



T/ \ nil /3+1 



(8) 



(9) 



(10) 



Here, Ai and Az are functions of temperature. 

In order to obtain an extended LIR EOS, we would limit parameters mi and ni to satisfy the relation- 
ship, mi =2 n\, and 

mi 13^ 2m, wi/3 = m. (11) 
Then, equation changes to the following form 

,2m+l 



RT 



-Ai 



+ Az 



(12) 



By using definition of compressibility, Z - PV I RT , equation (12) can be reformulated following the gen- 
eraUzed LIR (GLIR) EOS 

(Z-l) — =Bo + Bi y , (13) 



Bo 



AqVq 



A2V0 



(14) 



RT RT 

It can be seen that the LIR EOS in equation ([3) can be included in the GLIR EOS (13\ as a special case 
when m = 2. Since m = 1, equation (13) just reduces to the virial EOS. Although the parameter number of 
PSP EOS (D is the same as the three-parameter GLIR EOS (IS) , equation (13) with adjustable parameter m 
is more flexible and more accurate than equation l|4). 

Otherwise, we found in our calculations that the PM lllj and SSK (liKlll EOSs can be reformulated 
in the following forms: 



= D2 









Vol 


+ Co( 


-1 




Vol 








(vb 


\+Do 


{- 




[Vo 



(15) 



(16) 



We name this form as PMR and SSKR EOSs. Although the PMR EOS and SSKR EOSs are mathematically 
equivalent to the PM and SSK EOSs, they physically differ from each other. This is because all of equa- 
tions l[l), (|2) and equations (15) . (16) can be seen as Taylor expansion, but the expansion variable of equa- 
tions (D, (2) is (VqIV), and that of equations (15), (16) is (.V/Vq). At zero pressure, both values of (VqIV) 
and (V/Vq) are equal to 1. At high pressure, the values of (VIVq) are smaller than 1, the Taylor expansions 
in equations (15) . (16) are fast convergent. However, the values of (VqIV) are larger than 1 at high pres- 
sure, the Taylor expansions in equations (l), (l) are slowly convergent. Thus, the PMR and SSKR EOSs in 
equations (15) . (16) are more accurate than the original PM and SSK EOSs in equations (l), G). 
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3. Results and discussion 

Now we apply six EOSs to 28 metaUic solids, including GLIR |[13), PM (H, PMR (IS), SSK IH [ll, 
SSKR (16) and PSP yj] EOSs. All experimental data are taken from Kennedy and Keeler (1972) except 
forWl22J. 

Table 1. The experimental data of Vq (cm^/mol) and comparison of average relative errors (Ap%) of pres- 
sure for 30 metallic solids calculated from the GLIR, PM, PMR, SSK, SSKR, and PSP EOSs. 





Vo 


GLIR 

Ap% 


PM 

Ap% 


PMR 

Ap% 


SSK 

Ap% 


SSKR 

Ap% 


PSP 

Ap% 


Cu 


7.115 


0.54 


0.54 


0.50 


9.34 


4.65 


0.43 


Mo 


9.387 


0.80 


1.99 


1.45 


1.08 


1.05 


1.32 


Zn 


9.166 


0.30 


0.55 


0.39 


9.34 


5.46 


0.48 


Ag 


10.27 


0.38 


0.47 


0.41 


6.22 


4.10 


0.45 


Pt 


9.098 


0.70 


0.70 


0.70 


2.15 


1.69 


0.70 


Ti 


12.01 


0.68 


3.38 


2.03 


2.07 


1.07 


2.09 


Ta 


10.80 


0.66 


1.14 


0.86 


0.68 


0.64 


0.89 


Au 


10.22 


0.64 


0.64 


0.64 


2.44 


1.80 


0.64 


Pd 


8.896 


0.72 


0.72 


0.72 


2.21 


1.55 


0.72 


Zr 


14.02 


0.62 


7.02 


4.21 


3.70 


2.13 


3.51 


Cr 


7.231 


0.99 


1.00 


1.00 


1.35 


1.20 


1.00 


Co 


6.689 


0.65 


0.61 


0.61 


0.60 


0.60 


0.60 


Ni 


6.592 


0.60 


0.61 


0.61 


1.13 


0.98 


0.62 


Nb 


10.83 


1.71 


1.90 


2.00 


2.15 


2.20 


1.76 


Cd 


13.00 


0.24 


0.30 


0.29 


4.08 


2.92 


0.31 


Al 


10.00 


0.85 


0.76 


0.76 


0.87 


0.66 


0.49 


Th 


19.97 


0.36 


0.95 


0.69 


0.78 


0.53 


0.86 


V 


8.365 


0.39 


0.70 


0.58 


0.41 


0.39 


0.49 


In 


15.73 


0.57 


0.76 


0.62 


3.98 


2.91 


0.60 


Be 


4.890 


0.43 


0.64 


0.55 


0.61 


0.52 


0.51 


Pb 


18.27 


0.30 


0.31 


0.31 


2.86 


2.06 


2.09 


Sn 


16.32 


0.26 


0.31 


0.29 


2.50 


1.81 


0.29 


Mg 


14.00 


0.33 


0.59 


0.47 


0.27 


0.26 


0.66 


Ca 


26.13 


0.61 


5.67 


4.05 


2.91 


1.69 


4.98 


Tl 


17.23 


0.29 


0.29 


0.29 


1.38 


1.08 


0.28 


Na 


23.71 


0.48 


0.94 


0.68 


0.23 


0.21 


1.15 


K 


45.62 


0.39 


1.30 


1.20 


0.85 


0.46 


1.80 


Rb 


56.08 


0.43 


1.21 


1.15 


0.64 


0.36 


1.75 


mean error 


0.57 


1.29 


1.00 


2.39 


1.61 


1.12 



In table[T] we list the volume at zero pressure Vq, average fitting errors of pressure for the 28 solids. It 
can be seen that the GLIR (13) yields the smallest fitting errors for 20 solids, and for the other 8 solids the 
errors are also fairly small. The fitting precision for different solids is fairly stable for the GLIR EOS (13) . 
while instable for the other five EOSs. The largest errors among the 28 solids for the six EOSs are 1.71% 
of Nb, 7.02% of Zr, 4.21% of Zr, 9.34% of Zn, 5.46% of Zn, 4.98% of Ca, respectively. In the last line of the 
table, we list the total average error for the 28 solids. It can be seen that the GLIR EOS yields the best 
results with average error 0.57%; the PMR EOS yields second best results with average error 1.00%; the 
PSP EOS, PM EOS, SSKR EOS, and SSK EOS subsequently give worse results with average errors 1.12%, 
1.29%, 1.61% and 2.39%, respectively 

In tables[2]and[3l we list the fitted parameters for the six EOSs, table[2]shows that the values of m in the 
GLIR EOS (13) are smaller than 1 for 19 sohds, and slightly larger than 1 for 10 sohds. This imphes that the 
interactions in the metals are far softer than the LJ (12-6) potential, and are approximately approaching 
the LJ (6-3) potential for the 10 soUds, and even softer than the LJ (6-3) potential for other 20 solids. The 
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Table 2. Optimized values of coefficients for the GLIR, PM and PMR EOSs determined by fitting experi- 
mental compression data. The parameters for the GLIR EOS are dimensionless; and all parameters for 
PM and PMR EOSs are in GPa. 





GLIR 


PM 


PMR 




171 








^\ 


^2 








Cu 


0.906 


-449.53 


448.57 


-153.86 


167.29 


-13.44 


-13.24 


166.74 


-153.51 


Mo 


0.592 


-1731.45 


1732.80 


-395.31 


458.15 


-98.61 


-101.51 


465.36 


-363.72 


Zn 


1.199 


-188.20 


187.21 


-48.41 


34.78 


13.53 


14.11 


33.18 


-47.34 


Ag 


1.197 


-368.37 


367.35 


-81.38 


55.72 


25.60 


26.14 


54.37 


-80.54 


Pt 


1.031 


-1000.33 


999.46 


-272.52 


263.83 


8.69 


9.54 


261.87 


-271.40 


Ti 


0.416 


-1185.34 


1184.03 


-126.25 


164.13 


-37.40 


-40.20 


171.78 


-131.30 


Ta 


0.527 


-1655.67 


1654.63 


-218.00 


369.91 


-88.75 


-90.88 


375.05 


-284.06 


Au 


1.004 


-762.52 


761.51 


-184.79 


183.79 


0.70 


1.52 


181.85 


-183.36 


Pd 


1.031 


-683.16 


682.19 


-188.56 


181.27 


7.29 


7.11 


181.68 


188.80 


Zr 


0.197 


-2763.5 


2762.6 


-122.61 


166.63 


-43.34 


-48.00 


179.09 


-130.68 


Cr 


0.924 


-609.57 


608.36 


-209.26 


226.66 


-17.46 


-17.58 


226.93 


-209.41 


Co 


0.729 


-733.72 


732.70 


-257.16 


318.67 


-61.49 


-61.29 


318.20 


-256.89 


Ni 


0.893 


-565.39 


564.35 


-211.11 


234.21 


-23.09 


-22.88 


233.72 


-210.83 


Nb 


0.582 


-1282.5 


1281.3 


-250.03 


331.05 


-81.09 


-82.18 


333.60 


-251.51 


Cd 


1.223 


-216.51 


215.61 


-36.99 


23.21 


13.76 


13.84 


23.00 


-36.86 


Al 


0.719 


-452.88 


451.48 


-106.33 


132.17 


-26.01 


-25.01 


129.71 


-104.84 


Th 


0.613 


-701.54 


700.18 


-65.47 


81.33 


-15.71 


-16.41 


83.24 


-66.73 


V 


0.569 


-939.67 


938.80 


-223.73 


293.70 


-69.84 


-70.98 


296.38 


-225.29 


In 


1.058 


-240.83 


239.89 


-38.61 


36.62 


1.930 


2.430 


35.31 


-37.75 


Be 


0.477 


-500.16 


499.19 


-177.73 


239.41 


-61.58 


-62.52 


241.66 


-179.05 


Pb 


1.022 


-323.72 


322.61 


-43.61 


42.49 


1.090 


1.240 


42.12 


-44.37 


Sn 


1.118 


-256.62 


255.80 


-37.87 


32.02 


5.850 


5.970 


31.72 


-37.69 


Mg 


0.592 


-338.70 


337.60 


^4.32 


55.83 


-11.41 


-11.81 


56.91 


^5.03 


Ca 


0.076 


-2812.2 


2810.9 


-20.35 


27.60 


-6.810 


-8.040 


31.34 


-23.05 


Tl 


1.154 


-216.85 


215.66 


-28.85 


21.85 


6.990 


7.160 


21.43 


-28.60 


Na 


0.540 


-112.14 


111.09 


-6.740 


8.370 


-1.540 


-1.660 


8.740 


-7.030 


K 


0.419 


-147.58 


146.11 


-2.900 


3.670 


-0.670 


-0.750 


3.980 


-3.160 


Rb 


0.457 


-115.00 


113.29 


-1.780 


2.230 


-0.360 


2.440 


-1.970 


1.210 



table also shows that the parameter Bi in the GLIR EOS JT3) always takes on positive values, and this 
ensures a correct tendency as the volume tends to infinity. 

However, the values of C2 in the PM and PMR EOSs are negative for 18 and 25 solids, respectively. 
The values of D2 in the SSKR EOS, A2 in the PSP EOS are also negative for 2 and 18 solids, respectively. 
For these solids, the corresponding EOSs may exhibit a physically incorrect tendency as the volume tends 
to infinity. To compare, the GLIR EOS l[13) is not only the most precise one, but also is a unique EOS that 
does not exhibit a physically incorrect tendency among the six EOSs studied in this work. In figures [l] [H 
and[3l we plot the experimental compression data and the curves calculated using the GLIR, PMR, SSKR 
and PSP for 10 soUds, including Cu, Mo, Ag, Ti, Ta, Zr, Ni, Nb, Th and Be. These figures show that the 
calculated compression curves from the GLIR and SSKR EOSs are correct at high pressure for the 10 
solids, although for Zr, the parameter D2 in the SSKR EOS takes on a negative value. But the PMR EOS 
jdelds incorrect compression curves at high pressure for 7 of 10 solids, except for the solid Cu, Ag, and 
Ni. And the turn point is in the range VjVa =; (0.3 h- 0.5) for the 7 solids. Moreover, the PSP EOS also yields 
incorrect compression curves at high pressure for 9 of 10 solids, except for the solid Ag. And the turn 
point is about VIVa =; 0.3 for solids Cu and Ni; about VIVq =; 0.5 for other 7 solids. 

In these figures, we also plot the variation of relative errors of pressure with compression ratio VIVq. 
It can be seen from these figures that, for solids Cu and Ag, the oscillations of relative errors from the 
SSKR EOS are the most prominent, and are the same from other three EOSs; for solids Ti and Zr, the 
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150r 
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(f) 



Figure 1. (Color online) Comparison of compression curves of Cu (a), Mo (b), Ag (c), Ti (d), Mg (e), and Zr 
(f) calculated by using different equations with experimental data (o): solid line, PSP EOS; dashed line, 
SSKR EOS; dot line, PMR EOS; dot-dashed line, GLIR EOS. And comparison of percentage error of pressure 
calculated using different equations: +, PSP EOS; □, SSKR EOS; PMR EOS; O, GLIR EOS. 
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Table 3. Optimized values of coefficients for tlie SSK, SSKR and PSP EOSs determined by fitting experi- 
mental compression data. The parameters for the PSP EOS are dimensionless; and aU parameters for SSK 
and SSKR EOSs are in GPa. 





SSK 


SSKR 


PSP 




Do 


Di 


D2 


Do 


Di 


D2 


^0 


Al 




Cu 


293.73 


-686.40 


394.71 


381.53 


-650.37 


269.96 


427.41 


-419.20 


-9.020 


Mo 


156.70 


-578.16 


421.54 


421.11 


-577.10 


156.06 


1316.8 


-1209.8 


-107.85 


Zn 


260.70 


-528.62 


271.29 


256.94 


-489.17 


234.34 


192.99 


-206.99 


12.750 


Ag 


326.55 


-720.31 


395.34 


382.78 


-688.72 


307.02 


359.42 


-389.56 


28.99 


Pt 


466.76 


-1193.8 


727.49 


720.07 


-1176.7 


456.98 


996.17 


-1008.9 


11.81 


Ti 


8.260 


-121.04 


112.44 


114.63 


-127.03 


12.22 


585.21 


-536.06 


-49.24 


Ta 


72.99 


-344.84 


271.82 


272.46 


-346.38 


73.90 


1179.1 


-1066.7 


-113.03 


Au 


325.64 


-816.00 


490.93 


484.10 


-799.81 


316.14 


761.26 


-763.39 


1.120 


Pd 


344.40 


-865.63 


521.72 


513.51 


-846.47 


333.32 


682.39 


-690.84 


7.450 


Zr 


-40.14 


-22.08 


61.86 


64.60 


-29.40 


-35.39 


690.45 


-617.45 


-73.00 


Cr 


259.25 


-700.19 


441.11 


437.08 


-609.92 


253.56 


605.12 


-591.28 


-15.00 


Co 


161.67 


-516.19 


354.61 


354.11 


-515.05 


161.00 


663.78 


-618.34 


-46.25 


Ni 


241.20 


-661.39 


420.41 


417.51 


-654.68 


237.36 


553.90 


-537.02 


-17.83 


Nb 


70.990 


-313.58 


242.46 


241.83 
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(b) 

Figure 2. (Color online) Comparison of compression curves of Ni (a) and Na (b) calculated by using differ- 
ent equations with experimental data (o): solid line, PSP EOS; dashed line, SSKR EOS; dot line, PMR EOS; 
dot-dashed line, GLIR EOS. And comparison of percentage error of pressure calculated using different 
equations: -l-, PSP EOS; □, SSKR EOS; PMR EOS; O, GLIR EOS. 
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oscillations of relative errors from the PSP EOS and PMR EOS are more evident than the SSKR and GLIR 
EOSs; and for other solids, the oscillations from all four EOSs are equivalent with each other. It is notable 
that the relative errors from the GLIR EOS are most stable and fairly small for all 10 solids and for all 
compression ratio ranges. These results show that the GLIR EOS can be seen as the best one among six 
EOSs studied in this work. 




Figure 3. (Color online) Comparison of compression curves of Th (a) and Be (b) calculated by using differ- 
ent equations with experimental data (o): solid line, PSP EOS; dashed line, SSKR EOS; dot line, PMR EOS; 
dot-dashed line, GLIR EOS. And comparison of percentage error of pressure calculated using different 
equations: -l-, PSP EOS; □, SSKR EOS; PMR EOS; O, GLIR EOS. 



4. Conclusion 

In conclusion, we develop a three-parameter GLIR based on the GLJ potential and the approach of 
Parsafar and Mason 1 14] in developing the LIR EOS. Comparing with other five EOSs popular in literature, 
the precision of the GLIR EOS developed in this paper is superior to other EOSs. The GLIR EOS is capable 
of overcoming the problem existing in other EOSs where the pressure becomes negative at high enough 
pressure conditions. 
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Ysara/ibHeHe pibhuhha craHy, sacTOCOBHe f\o Mera/iiB 

r. Cai^, AxJ. Ca^P^, B.flx. M^, flx. TaHi^ 

Katjieflpa npuK/iaflHOi (fiiswKH, KwTaticbKWM yHisepcuTeT e/ieKipoHiKH ra TexHO/ioriti, MeHfly 610054, KHP 
/la6opaTopifl c|)i3HKH yflapHoTxBH/ii i fleTOHai^ii, niBfleHHO-SaxiflHuPi iHCTwryT c|)i3HKH n/inHlB, 
MiaHSH 621900, KHP 

BanponoHoeaHO TpunapaMerpuHHe pisH^iHH^i craHy 5e3 c|)i3H4H0 HeKopeKTHi/ix oci4n/iflL4iR, ll(0 6a3yeTbCfl Ha 
y3ara;ibHeH0My noTeHi^ia/ii /leHHapfla-flxoHca (GLJ) i niflxofli napcac|)apa i MeticoHa [Parsafar G.A., Mason E.A., 
J. Phys. Chem., 1994, 49, 3049] flo BUBefleHHfl piBH^HHA craHy 3 pery/iapHiCTK) yilHiwHOi isorepMH (LIR). Sanpo- 
noHOBane ysara/ibHeHe pIbhahha craHy Moxe BK;iK)4aTH b ce6e LIR piBH^inna CTany nacTKOBuPi BunaflOK. 
Tpi/inapaMeTpn4Ha ysara/ibHena pery/iapHiCTb /imiCiHoT i30TepMW [Parsafar G.A., Mason E.A., Phys. Rev. B, 1994, 
49, 3049] (PM), [Shanker J., Singh B., Kushwa S.S., Physica B, 1997, 229, 419] (SSK), [Parsafar G.A., Spohr H.V., 
Patey G.N.,J. Phys. Chem. B, 2009, 113, 1 1980] (PSP) i nepec|)opMy;ibOBaHi PM SSK piBHflHHfl ciaHy e 3acT0C0Ba- 
Hi flo 30 Mera/iiHHWx TBepflux ri/i y LuwpoKiti o6;iacTi rwcKy. noKa3aHO, L140 PM, PMR i PSP piBHflHHa craHy fl/ia 
6i;ibLiJ0CTi TBepflwx tI/i ra SSK i SSKR piBH5iHH5i CTany fl/ia fleKi/ibKOx iBepflux ri/i Maroib c|)i3M4H0 HeKopeKThi 
noBopoTHi TOHKn, i TMCK ciac HeraiwBHWM npn flocnTb BucoKOMy TMCKy. y3ara;ibHeHe pIbhahha craHy e 3fla- 
THWM He Ti/ibKW HOflO/iaTM npo6;ieMy, icnyiOHy b Ihiumx n'arn pIbhahhax crany, fle tuck crae HerarnBHUM npn 
BucoKOMy TWCKy, a/ie raKOx flae Kpaini pesy/ibrarw, nix miui pIbhahha crany. 

KniOMOBi cyiOBa: Tpi/inapaMeTpi/mne pibhrhhr crany, MeTani^ni TBepfli ri/ia, BucoKud ti/ick, cpisn'^HO 

HeKOpeKTHi OCU,l/inRU,li 
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